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Abstract 

In this paper we examine the spatio-temporal dynamics of two nonhnearly 
coupled wave triplets sharing two common modes. Our basic findings are 
the following. When spatial dependence is absent, the homogeneous manifold 
so obtained can be chaotic or regular. If chaotic, it drives energy diffusion 
from long to small wavelengths as soon as inhomogeneous perturbations are 
added to the system. If regular, one may yet have two distinct cases: (i) 
energy diffusion is again present if the inhomogeneous modes are linearly 
unstable and triplets are effectively coupled; (ii) energy diffusion is absent if 
the inhomogeneous modes are linearly stable or the triplets are uncoupled. 
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I. INTRODUCTION 



Wave interaction in continuous media is frequently described in terms of amplitude equa- 
tions. The procedure used to obtain these amplitude equations is standard. Once one iden- 
tifies a set of high-frequency modes interacting in a resonant fashion, one applies multiple 
scale formalisms to derive equations governing the slow space-time evolution of these modes; 
the amplitudes equations. 

One of the most investigated types of resonant interaction that can be effectively de- 
scribed by amplitude equations corresponds to the triplet interaction. In that case, three 
high-frequency modes undergo a resonant interaction which, in the absence of external drives 
and dissipation, can be shown to produce integrable space-time dependent amplitude equa- 
tions tl. 

Another case of resonant interaction that has been attracting some renewed attention is 
the one involving two coupled triplets sharing two common modes This type of in- 

teraction has been shown to produce amplitude equations allowing for temporal chaos . 
To be more specific, if one removes all spatial dependence from the governing equations, 
then the solutions developing on this homogeneous manifold can be chaotic if a number 
of conditions are fulfilled. Not much, however, has been said about solutions departing 
from the homogeneous manifold. One relevant question here would be on the stability of 
the homogeneous manifold against inhomogeneous perturbations. If, for instance, a cer- 
tain homogeneous solution - not necessarily a static equilibrium - is perturbed with some 
small inhomogeneity, would this inhomogeneity grow in time? If the answer is positive we 
would have an unstable homogeneous manifold. Otherwise the manifold would be stable or 
marginally stable. One could actually split the stability analysis into the analysis of ho- 
mogeneous manifolds supporting regular dynamics and homogeneous manifolds supporting 
chaotic dynamics, and the respective results could be significantly different. In fact, as we 
have observed and shall discuss in detail later, chaotic manifolds are intrinsically unstable. 
In other words, while regular manifolds may be unstable against some types of perturbation 
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and stable against others, chaotic manifolds always tend to be unstable no matter the type 
or strength of the perturbation. What happens in the latter case is that inhomogeneities 
always tend to grow as a result of the stochastic drive provided by the chaotic orbits evolving 
on the homogeneous manifold. Here one cannot use linear techniques to investigate stability 
since the problem is intrinsically nonlinear; one should make use of stochastic models 
instead. On the other hand, if homogeneous chaos is absent the stochastic drive is absent 
as well. Then stability can be estimated by linear analysis which in fact shows that only 
modes falling within the linear instability band are the ones to grow. 

Based on these comments one can already see that at least in the case of chaotic homoge- 
neous manifolds, the instability is unlikely to be saturated by any nonlinear effect involving 
a small number of modes, like soliton formation for instance. As energy is continuously 
fed into the inhomogeneities, modes with smaller and smaller length scales are eventually 
excited. Interestingly, we shall also see that in the case of regular homogeneous manifolds 
coupled to linearly unstable modes, the behavior is similar; the initial instability is not satu- 
rated by nonlinear effects involving a finite number of modes and once more energy spreads 
toward smaller and smaller length scales. Energy flow is arrested only when the triplets are 
uncoupled, or else when the homogeneous manifold is regular and no inhomogeneous mode 
falls within the linear instability band. 

In view of the previous facts, the issues to be discussed here are related to topics on 



turbulent motion or on the Fermi-Pasta-Ulam problem |1T0|-[12|. Indeed, our homogeneous 
manifold could be seen as representing fluctuations with very long wavelengths. Therefore 
the stability problem that we wish to pose here could be more precisely stated as follows: 
when is our system unstable and under which conditions can energy flow from long to short 
wavelengths? Contrarily to integrable cases where no flow is observed, we shall see here that 
energy transfer or energy redistribution does occur under the conditions mentioned above: 
(i) when the homogeneous manifold is chaotic; or else (ii) even with a regular homogeneous 
manifold, when the two triplets are effectively coupled and the inhomogeneous perturbations 
fall within the instability band. 



The paper is organized as follows: in §11 we introduce the basic formalism, governing 
equations, and numerical methodology; in §111 we analyze the dynamics on the homogeneous 
manifold with help of Poincare maps, in §IV we perform the relevant analyses and simulations 
for the full spatio-temporal problem, and in §V we conclude the work. 

II. THE MODEL 

We are interested in describing the mutual interaction of two wave triplets sharing two 
common modes and constituted, therefore, by four wave modes in all. We assume that the 
following resonant conditions are fulfilled: 

Kl = K4 - K2, (1) 

where cu^^ and k^, {j = 1,2,3,4}, are respectively the fast frequencies and fast wavevectors 
of the interacting waves. 

Note that we allow for small frequency mismatches, an usual effect in wave-wave in- 
teraction. Frequency mismatch occurs because even for perfectly matched wavevectors, the 
respective frequencies obtained from the linear dispersion relations may not be quite likewise 
matched. In laser accelerators and in general laser-plasma interactions, frequency mismatch 
can even enhance the linear growth rate in several situations [f7i[T3|. Therefore, mismatched 
modes can be stronger and of greater importance in the dynamics. Moreover, as we shall see 
here, the integrability properties of the nonlinear regime of the interaction depend critically 
on the existence of mismatch. Wavevector mismatches could be also incorporated into the 
theory. However, as frequency and wavevector mismatches are formally equivalent we focus 
attention on the former. 
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This type of kinematics where the interacting triplets share two common modes, 1 and 
2, has been shown to be of relevance in nonlinear interactions involving electromagnetic, 
Langmuir and ion-acoustic waves propagating in unmagnetized plasmas [§,0, in magneto- 
hydrodynamics interactions involving Alfven and ion-acoustic waves 0, and in plasma-beam 
interactions in the presence of negative energy waves Following the model revisited in 
a series of recent papers the dimensionless amplitude equations governing one dimen- 
sional, spatio-temporal, slow modulational dynamics can be cast into the form: 

+ = Mx, t)A,{x, t) - r A*{x, t)A,{x, t), (2) 

+ ^9. = -M^, t)A;{x, t) - r Alix, t)A,{x, t), (3) 

^ +vg, g^' ^ = i6sAs{x,t)-Ai{x,t)A;{x,t), (4) 

dAAx,t) dAJx,t) .r^,^ . / X . / X /^S 

' +Vg, ' = i6iA^{x, t)+r Ai{x, t)A2{x, t). (5) 

Aj, {j = 1, 2, 3, 4}, are the complex amplitudes of the four fields, ^3^4 = ui ^ 002 — 4 are 
independent frequency mismatches corresponding to fields ^3 and A4 (i.e., one can always 
take 61 = 62 = 0, as we actually did), r is a variable strength factor measuring the intensity 
of the triplet-triplet coupling ^j^, and Vg^ are the respective group velocities along the 
spatial modulation that we take as one dimensional (the x axis) in this work. Time and 
space derivatives are first order as a result of our multiple time and space scales. 

The set of governing equations (|D-(|5D can be derived from a continuous Hamiltonian. 
Indeed, it is possible to see that the following relations hold 

dAj{x,t) _ SH dAj{x,t)* _ SH 



dt 6 A*' dt 6Aj 
where we introduce the functional derivative 

5 d d d 



(6) 



SAj dAj dx 5(^) 



(7) 



(the same definition holds if A is replaced with A*), and where the Hamiltonian is to be 
written in the form 



H = J dx[-A,A*Al + AIA2A3 - r{AlA'^A^ - A,A2Al)+ 

dA^ 
dx 



^S^\A,\^ + ^5M^\'-E^^.A^^]■ («) 



The Hamiltonian does not depend explicitly on time; therefore it is a time conserved quantity. 
In addition to the Hamiltonian, the following quantities are also conserved: 

Ci = J rfx[|A2|'-|A3|2 + |A4|'], (9) 

C2 = Jdx + + , (10) 

and they suggest that we can look at the whole process as a decay interaction with Ai as the 
decaying pump; from this perspective, A2 is an idler wave whereas ^3 is a Stokes mode and 
A4 an anti-Stokes mode. In the case of nonlinear interactions in unmagnetized plasmas for 
instance, wave Ai is a transverse electromagnetic wave, wave A2 is ion-acoustic, and waves 
A3 and are Stokes and anti-Stokes Langmuir modes 

All those quantities will be useful in checking out the accuracy of our integration methods 
which we outline now. The basic integration method is pseudo-spectral and the four fields 
Aj , {j = 1, 2, 3, 4}, are Fourier analyzed according to 

N/2 

Aj= E «..We^"'^ (11) 

n=-N/2+l 

For further purposes note that while j denotes the wave type, n denotes the mode number 
or harmonic number. The harmonic or modal amplitudes aj^ are integrated in time with 
a predictor-corrector algorithm. We use = 64, 128, 256 modes removing half of them 
to cure aliasing. We denote the basic slow wavevector by k and point out that due to 
the structure of the equations, variations in k can be absorbed in variations of the group 
velocities or vice-versa. Fluctuations of the conserved quantities, including energy, are not 
larger than one part in 10^, and variations of the tolerance factors of the integrator do not 
alter the final outcome of the runs. 
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III. DYNAMICS ON THE HOMOGENEOUS MANIFOLD 



As a first step let us analyze the behavior of the dynamics on the homogeneous manifold. 
The motivation and relevance of a consistent analysis of the homogeneous dynamics has been 
commented in the Introduction and has roots in the stochastic pump model that we discuss 
next. According to the stochastic pump model 0, as soon as strongly chaotic degrees-of- 
freedom of a dynamical system are coupled to additional degrees-of-freedom, the chaotic 
subsystem acts like a stochastic pump delivering a net amount of energy to its environment. 
This net energy transfer is diffusive and takes place at all because of the very random 
nature of the coupling. In a certain sense, the dynamics is equivalent to the coupling of 
two thermodynamical systems, one cold and the other hot. As a result of the random 
coupling, a net energy transfer finally equalizes the temperatures of both systems. In our 
case, the hot system would be the chaotic homogeneous manifold and the cold one the 
subset of the inhomogeneous modes. Furthermore, and precisely due to diffusion, energy 
would be redistributed over the accessible phase-space of the inhomogeneous modes exciting 
progressively harmonics with finer and finer length scales. 

Homogeneous orbits have been recently investigated in a paper by Pakter, Lopes, & 
Viana Let us review the basic results and take the opportunity to introduce a convenient 
notation to be used here. To investigate the types of orbits on the homogeneous manifold, 
let us thus assume that we are working with fields of the form 

A,{x,t) = ajSt), (12) 

which do not depend on the spatial coordinate. Then, introducing real amplitudes and 
phases through 

ajo = VPje'^^, {j = 1,2,3,4}, (13) 

one obtains canonical conjugate equations for these quantities, with a reduced governing 
Hamiltonian h given in the form 
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h = 2y/pip2P3 sin(0i - 02 - 03) - 2r^pip2P4 sin(0i + 02 - (pi) - ^sPa - S^p^- (14) 
This type of Hamiltonian structure has been largely studied in the context of nonlinear 



waves 



p!4| . What we learn is that due to the particular phase combinations of the sine 
functions, a number of relatively simple canonical transformations can be applied to simplify 
the problem. The first transformation introduces a new phase angle in the form = 
01 — 02 — 03 and preserves the other phases, 0^- = 0j, {2,3,4}. Rewriting conveniently the 
argument of the second sine function of expression (p!^), the form of the second canonical 



transformation becomes clear; it shall be constructed such as to introduce a new phase angle 
04 = 04 — 202 — 03, preserving all the others 0^' = 0^-, {j = 1, 2, 3}. As a result of the pair 
of transformations described above, one finally arrives at 

P2 = p2 - Pi - 2P4, (15) 
P3 = P'3' - P'l' - pI (16) 

and 



h = 2Jpi{p2 - pi- 2p4){p3 - pi - P4) sin( 



2rJpiPi{p2 -pi- 2p4) sin(0i - 04) - 5s{p3 - pi- p^) - S4P4, (17) 



where the double primes have been dropped in the Hamiltonian (^7^. Since the Hamiltonian 
no longer depends on coordinates 02 and 03, the associated momenta p2 and p3 are constants 
of motion with their relation to the original variables given by relations (|15]) and (|16|). 



Furthermore, the Hamiltonian h, Eq. (p!?!), is a conserved quantity itself since it does 
not depend explicitly on time. We now focus attention on conditions allowing for decay 
processes where Ai is to be considered a pump wave. In the decay case one requires Ai ^ 
and A2 = A^ = A4 = 0, which implies, in view of conditions (|12|)-(|l7|), that our energy 
surface of interest here is actually indexed hj h = 0. 

The flow governed by h is two-degrees-of-freedom and thus likely to be nonintegrable . 
The relevant governing equations obtained from Hamiltonian (|T7|) are: 
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^ = -2y'pi(p2 - pi- 2pi){p3- pi - Pi) cos 01 



+ 2rJpipi{p2 - pi- 2pi) cos( 



dpi 
dt 



-2rJ pipi{p2 - pi- 2pi) cos( 



(18) 
(19) 



dt 



\jpl{.p2 - Pi- 2P4)(P3 - Pi - P4) 



sm 01 X 



P2P3 - 2pip2 - P2P4 - 2piP3 + 3pi + 6pip4 - 2p4p3 + 2p^ 



and 



r sm(0i - 04) 



PiP4(p2 - Pi - 2p4) 



P4P2 - 2piP4 - 2p^ + ^3, 



(20) 



dt 



Pl(P2 - Pi - 2P4)(P3 - Pi - P4) 



: sm 01 X 



-2pip3 + 3pi + 4pip4 - P1P2 



r sm(0i — 04 1 



P1P2 - Pi - 4piP4 



+ ^3 - 



(21) 



'PlPi{P2 - Pi- 2p4 

However there is at least one pair of limiting cases where the dynamics becomes com 



pletely integrable. One of them occurs when ^3 = 54 = 0. In this situation Romeiras |15 



has identified a new constant of motion that reduces the system to one-degree-of-freedom. 
The other integrable limit occurs when r — 0; in this situation mode decouples from the 
multiplet and one recovers the pure three wave decay which, again, is one-degree-of-freedom 
and hence completely integrable. 

To understand and visualize the properties and characteristics of the purely homogeneous 
interaction, let us fix r = 1.0 - a typical condition of interactions processing in unmagnetized 
plasmas or magnetohydrodynamics environments, for instance - and ^4 = 0. Let us then 
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proceed to investigate the system behavior as S3 is allowed to vary. We make use of Poincare 
plots, where we record the values of pi and 0i, as determined by the governing reduced 
Hamiltonian (0), each time dp^/dt = with d^pi/dt^ > 0. We consider weak nonlinearities 
and launch a number of initial conditions, in fact 10, under the constraint h = 0. In all the 
conditions we take 0i = 71/2, 04 = 0; pi is uniformly distributed within a weakly nonlinear 
range < pi < pmax < 1 and p4 is calculated from h = 0. Finally, to be consistent with 
typical decay conditions, we take P2 = Ps = Pmax such that whenever pi = p^ax and p4 = 0, 
A2 = A3 = 0. 

As mentioned, in this paper we investigate weakly nonlinear regimes; all the simulations 
are performed with Pmax = 0.1. Now the first plot shown in Fig. (|l]a) is produced for 
parameters lying in an essentially regular regime where r = 1, and ^3 = 0.0001. As predicted 
by Romeiras [|15| , for such a small ^3 the phase-space is essentially regular displaying a series 
of closed curves that in fact reveal the existence of numerous invariant KAM (Kolmogorov- 
Arnold-Moser) tori 0. Then, as we increase the value of the mismatch up to ^3 = 0.1, Fig. 
(|I|b) indicates that the KAM tori become substantially eroded; chaos is strong in this case. 
As we couple the homogeneous manifold to inhomogeneous perturbations, diffusion towards 
small wavelengths is expected to occur in the latter situation. 

We have therefore obtained some initial information on the behavior of the homogeneous 
manifold. But the real question yet to be answered is how that kind of behavior on this 
manifold affects the inhomogeneous dynamics. We shall address this issue next. 



IV. FULL SPATIO-TEMPORAL DYNAMICS 

The typical initial conditions to be used in this paper place the homogeneous orbits on 
the outermost curve of Fig. (|l]a), close to a separatrix orbit containing two unstable fixed 
points located at the maximum pi - in this case pi = 0.1 - and (pi = 0, tt. These fixed points 
represent the exact homogeneous equilibria of our system. In regular cases the orbit remains 
close to the separatrix and since the separatrix orbit spends an infinitely large amount of 
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time close to the fixed points, we estimate linear stability simply by replacing the actual 
orbit with the fixed points. The estimates shall be shown to be accurate enough when 
compared with the simulations. 

Let us thus suppress for a little while the dynamics on the homogeneous manifold and 
assume that we have an equilibrium situation characterized by A^'^ = af'^ = ^ 7^ is a 
complex constant), and afj = , {j = 2, 3, 4}, for any other n. As mentioned before, this is 
the convenient condition if one wishes to study the stability of a strong pump - mode Ai - 
feeding vanishingly small modes - modes Aj, {j = 2,3,4}. The stability of the pump can 
then be readily investigated if one adds small perturbations A^j^^ = af^e^^^, a^^^ ~ ^-tnt 
the system. Indeed, if {af^l <^ l^i^'*!, a dispersion relation = ^l{k) can be obtained from 

P3P2-Pr = 0, (22) 

where 

P3 = + Vg^k - 53), 

Pn=\Ani-A 
-T 4 

P4 = + Vg^k - 54). 

The dispersion relation (^2]) can be analyzed numerically but here we abbreviatedly discuss 
some of its relevant results. For a given A, there exists an instability (complex f2's) band 
extending from A; = up to a certain k^ax- The larger the value of A, the larger k^ax, but 
otherwise the instability band is not strongly dependent on ^3, ^4, and r. The instability for 
= 0, in particular, is the one giving rise to the homogeneous dynamics that was the subject 
of study in section II and in previous works As a relevant information obtained there, 
we recall that the ensuing homogeneous dynamics can be chaotic or regular depending on 
the initial conditions and control parameters. 

We can now collect the results obtained so far to assert that in general we have four 
global types of situations to be investigated. These four situations result from the four 
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distinct associations involving a regular or chaotic homogeneous manifold, with inhomo- 
geneous perturbations launched inside or outside the linear instability band. Our interest 
is to discover in which of the four situations can one detect energy diffusion from long to 
small wavelengths. Although all four cases are of interest, in this paper we shall focus on 
the two independent situations that provide a clear view on the individual role of the basic 
nonlinear effects driving diffusion: chaotic homogeneous manifolds associated with linearly 
stable inhomogeneous perturbations, and regular homogeneous manifolds associated with 
linearly unstable inhomogeneous perturbations. The behavior of the two remaining cases 
may be then obtained by an extension of the behavior of these basic ones. Our classification 
outlined above may sound artificial since the homogeneous and inhomogeneous dynamics 
are in fact coupled and mutually interactive; in other words they should not be considered 
apart from each other in principle. However, the classification adopted explains well enough 
what we are about to see: in the case of chaotic homogeneous manifolds coupled to stable 
modes, the diffusion is so slow that at any instant one can consider the chaotic dynamics as 
adiabatically delivering energy to the inhomogeneous modes; in the case of regular homo- 
geneous manifolds coupled to unstable modes, similarly, a possible regular dynamics on the 
homogeneous manifold does not alter significantly the results of linear stability calculations. 

To monitor instabilities and energy transfer, we make use of a spectral average |16|JT7[| 
that enables to estimate the number of active modes in the system. We denote this quantity 
by \/ < N"^ > and define it as 



average is taken over mode number {n) and mode type (j), and weighted by the square 
of the mode amplitude for a fixed initial condition. One could also average over several 
similar initial conditions but our results on energy transfer remain the same as long as the 
initial conditions are all simultaneously stable or all simultaneously unstable. v<iV^> is 
expected to grow in time in diffusive cases where more and more modes become involved in 




From its definition one sees that V < N"^ > is the square root of the averaged n^. The 



12 



the dynamics. In the absence of energy transfer, ^ < > remains hmited by the number 
of hnearly unstable modes; in the case of stable modes only, v<iV^> 0. 

In our simulations we take ^4 = 0. Our results on energy transfer are qualitatively inde- 
pendent of a precise choice of group velocities. Therefore we consider Vg^ = 0, {j = 1,3,4} 
for simplicity, since variations of Vg^ alone provides an easy way to control the width of the 
band of linear instability. We point out that a precise choice of velocities would be essential 
either in the case of a pure triplet interaction, or else if one were performing a detailed study 
of how the dynamics of the four-wave system evolves in time during intermediary stages, as 
the asymptotic state is approached. Indeed, depending on the ordering of group velocities 
solitons can actually saturate the pure triplet interaction or perhaps can serve as metastable 
intermediary states existing for a finite amount of time when a fourth wave destroying inte- 
grability is added. However, once the four-wave system has a chaotic homogeneous manifold 
or becomes linearly unstable as discussed earlier, energy transfer progresses regardless the 
presence or absence of solitons in these intermediary states. Note that when the coupling 
parameter r is sufficiently large, the dynamics of the four-wave system does not resemble 
that of the triplet interaction. This is why we adopt the simple ordering discussed above; 
other orderings not presented in the present work have been tested as well. 

In any case, such a choice is of physical significance for the interaction of electromagnetic, 
ion-acoustic and Langmuir waves in plasmas with small enough Debye length. In that case 
one shows that the only relevant group velocity is the one corresponding to the ion-acoustic 
mode, the A2 wave 0. The choice of the numerical value for Vg^ and k is arbitrary to a large 
extent. This is due to the structure of the governing equations (0)-®, which are invariant 
under the rescalings 

Vg^ ^ avg^, X ^ ax (24) 

and 

A, PA„ t i, X ^ ^, 5, ^ (36„ {j = 1, 2, 3, 4}, (25) 
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a and /3 being independent scale factors. This fact effectively creates some additional free- 
dom in choosing the simulation parameters. Here we shall consider /c = 1.0 in all cases. 



A. The ^3 7^ case. 

Let us consider as a first instance the case of a chaotic homogeneous manifold combined 
with linearly stable perturbations. To simulate this situation, we use the same parameters 
and techniques employed to construct Fig. (Pd). In the homogeneous manifold the initial 
condition corresponds to the outermost orbit present in the figure. Since we had launched 
10 equally spaced orbits, a rapid calculation therefore shows that the value of p\{t = 0) 
to be used in the simulation reads pi(t = 0) = 0.09. Fig. (|l|b) also shows that we take 
01 (t = 0) = tt/2 for the initial phase. In addition we recall that 04 = and h = 0. Note 



that the relations between p's, 0's and aj^ are the same as previously defined in Eq. (|13|). 
For the inhomogeneous perturbation, besides k = 1.0, we take Vg^ = 1.0 as the relevant 
group velocity. The perturbing field itself is written in the form i?e[aij = |ai(j| x 10~^, 
/m[aij = |aig| x 10^'^ where i?e[«] and Jm[«] respectively denote real and imaginary parts 
of the various fields. For this particular choice of parameters and initial conditions, the 
inhomogeneous perturbations alone would be linearly stable, since here we have kmax ~ 
0.4 < k = 1.0. On the other hand the low- dimensional phase-space is chaotic, and because 
of this fact, some kind of energy transfer should be expected after all. The initial results 
of our simulations can be found in Fig. (^). It is clearly seen that in spite of the linear 



stability of the system, transfer does indeed occur as v<iV^> continuously grow until it 
reaches saturation due to the finite number of modes used in the numerical code. What 
happens here can be understood in terms of the stochastic pump model mentioned before, 
and discussed, for instance, in the book of Lichtenberg & Lieberman 0. The strong chaotic 
dynamics developing at = acts like a random force, diffusively driving all the remaining 
modes of the system. It does not matter if these remaining modes are linearly stable. What 
really matters is that as a result of the stochastic pump the inhomogeneous modes start to 



14 



be driven similarly to what would happen if they were in contact with a thermal reservoir. 
Simulations with other values of ^3 are also shown in Fig. (0). It is seen that as one lowers 6^, 
redistribution weakens. In the limit of very small values of the mismatch, one recovers the 
situation depicted in Fig. (|l|a) where the homogeneous manifold is regular. Then diffusion 
apparently ceases or becomes very slow, what leads us to draw a most expected conclusion 
on the features of the system: when the homogeneous manifold is regular and no harmonic 
mode falls within the linear instability band, diffusion is negligible. 

Figure (|^a) shows that the larger the number of modes used in the simulations, the 
larger the saturated value of V < N'^ >. In other words there is no nonlinear saturating 
effect that involves only a small and fixed number of modes. In addition to that, the various 
numbers of modes used in the figure, = 64, 128, and 256, indicate that diffusion is 
not an effect connected with the finite discretization of the numerical method. Diffusion 
is present no matter the value of A^. One can see this because for earlier times of the 
nonlinear interaction where a few modes are actually excited, the behavior of all curves are 
virtually the same. A close inspection of the final spectral-energy distribution reveals that 
our asymptotic states are not quite fully equipartitioned in energy. Although the spectral 
distribution for large mode numbers is relatively flat, the energetic content in the spectral 
region of small mode numbers is slightly larger. This can be seen with help of Fig. (^d) 
where we plot the energy of the free-modes. En =<^ {6j + nk) |aj„p doubly averaged over 
time and over the wave index {j = 1,2,3,4}. This discrepancy between full equipartition 
and simulations could have two contrasting sources that should be investigated: either the 
wave amplitudes we work with are already so high that the interaction term is no longer a 
small perturbation, or for the parameters used here integrable features are still appreciable. 
It is perhaps convenient to mention that the convergence of the spectral distribution has 
been verified by extending and reducing the averaging time interval in various ways; for 
sufficiently long intervals the results are the same. 

To complement the information obtained with help of Figs. and (^, in Fig. 
we display the space-time history of |74i(a;, t)| corresponding to the ^3 = 0.1 curve of Fig. 
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(0). It may be seen that for earlier times of the nonhnear interaction, one has a purely 
time dependent solution that is completely homogeneous in space. However, as time is let 
to evolve, strong spatial oscillations adds to the previous erratic but pure time dependence. 
The upper time limit of the simulation is t = 20000 which is smaller than the asymptotic 
time where the system reaches saturation. Nevertheless, it is clear from the figure that 
inhomogeneities are already strongly excited by then. 

B. The 63 = case. 

Now let us consider the other individually relevant situation in which we combine a 
regular manifold with a linearly unstable perturbation. The relevant graphical information 
can be found in Fig. (^). In this case we work with = so as to guarantee that the 
homogeneous manifold is in fact regular. The choice of initial conditions follows the pattern 
used to construct the previous figures. We launch 10 initial conditions on the homogeneous 
manifold equally spaced between pi = and pmax = 0.1, and select the initial condition for 
the full simulations as the one corresponding to the outermost orbit with (pi = tt/2. We 
choose = 0.06 such as to have at least some linearly unstable modes. In fact, within 
the range of variation of r in the figure, use of Eq. (^2|) shows that the number of unstable 
modes is approximately 5 for all cases. 

One could think that as chaotic activity is absent from the homogeneous manifold, no 
diffusion would be observed now. This proves to be untrue. Indeed, Fig. (1) shows that when 
r 7^ 0, energy diffusion is again present, even though homogeneous chaos is not (remember 
that 63 = 0). We had seen before that for coupled triplets - that is if r 7^ - purely 
temporal chaos alone is capable to drive energy redistribution, even though inhomogeneous 
perturbations are linearly stable. Now what we are observing is a complementary effect. 
Namely, unstable inhomogeneous modes alone are capable to drive redistribution if r 7^ 
even though the homogeneous manifold is integrable (^3 = 0). In this case the active modes 
cannot be excited with the stochastic pump simply because there is no stochastic pump 
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available. Instead they are excited by the linear instability of the system, and as soon as 
they are excited diffusion starts to operate again. This feature leads us to reconsider the first 
case analyzed here. In fact, in this first case where we have a chaotic homogeneous manifold 
combined with stable spatial modes diffusion can be understood in terms of the stochastic 
pump model only for earlier times of the nonlinear interaction, when the inhomogeneity is 
relatively small. Only then are the effects of the spatial inhomogeneities small, with diffusion 
totally dictated by the homogeneous chaos. As soon as the inhomogeneous fluctuations grow 
larger, both nonintegrable effects become of the same order of magnitude. At this stage the 
overall dynamics should become similar to that case where one combines a chaotic manifold 
with unstable inhomogeneous modes, a case not individually explored in this work. 



V. FINAL REMARKS 

A general comment is that in the wide sense investigated here, the global stability of the 
homogeneous manifold cannot be simply understood in terms of a linear stability analysis. 
In fact, we have shown explicitly that the linear analysis fails when the homogeneous man- 
ifold is chaotic. To emphasize our point we have analyzed a particular setting where the 
inhomogeneous modes that we have launched are linearly stable. Even so, inhomogeneous 
perturbations tend to grow as time evolves. What happens in this case is that at least 
for earlier times of the nonlinear interaction, the chaotic homogeneous manifold acts like a 
stochastic pump delivering energy in a random fashion to these inhomogeneous modes. One 
calls that situation the stochastic pump because the role of the homogeneous manifold is 
comparable to the role of a thermal reservoir diffusively delivering energy to its vicinities 

in- 

When homogeneous chaos is absent, two types of behaviors could be identified for r 7^ 
0. If no inhomogeneous perturbation falls within the band of linear instability, then the 
system does not develop any sort of configuration with spatial dependence and simply keeps 
supporting a homogeneous periodic dynamics. On the other hand, when at least one of the 
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inhomogeneous modes does fall within the instability band one has again diffusion towards 
small spatial scales. This shows that the nonintegrable features of our system are not 
exclusively related to the temporal dynamics developing on the homogeneous manifold. In 
the homogeneous manifold, chaos vanishes as 63 — > 0, but for dynamical states involving 
unstable modes with k ^ 0, energy redistribution is present irrespectively of the value of ^3. 
In other words, even if ^3 = in these situations, random activity is still present. 

The spatio-temporal dynamics of one isolated wave triplet has been shown to be inte- 
grable via inverse-scattering Our results indicate that as one couples a fourth wave 

such as to form two wave triplets sharing two common modes, integrability breaks down. 
What seems to be happening here is that as one couples the fourth mode, the regular so- 
lutions of the isolated triplet undergo a transition to chaos. Then, the resulting chaotic 
dynamics enhances radiative effects and causes energy transfer to modes with larger and 
larger wavevectors; this and the following are issues under current investigation. The regular 
asymptotic dynamics of the pure triplet is somewhat related to the intermediary dynamical 
stages of the four-wave system. Consider, for instance, an ordering of group velocities for 
which solitons saturate the pure triplet interaction, as commented earlier. When the fourth- 
wave is coupled to the unperturbed triplet, integrability is destroyed and solitons survive 
only for finite amounts of time, ultimately releasing their energy into small scale fluctua- 
tions. In other words, regular structures that are the final asymptotic states of the triplet 
interaction, are present only during intermediary stages of the full four-wave system. In 
addition, although the ordering of the group velocity is relevant during intermediary stages 
of the full interaction, we see that it is not expected to drastically affect the asymptotic 
dynamics and the overall process of energy transfer: eventually modes with short spatial 
scales are excited and intermediary regular structures are no longer present in the system. 

Another important point to be considered is that the general technique of multiple time 
scales leading to generic amplitude equations breaks down when modes with large enough 
wavevectors become involved in the dynamics. Then, dispersive and dissipative effects should 
be included since these terms are always related to small wavelengths - dispersion because it 
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appears in the form of derivatives terms like d'^Aj/dx^ such that the smaller the length scale, 
the larger the derivative, and dissipation because it occurs when the waves strongly interact 
with small scale structures of the global system, like for example high-velocity particles in 
plasmas. While we are currently considering these points, we emphasize that our present 
model works well for earlier times of energy transfer when modes with small wavelengths 
are still of very small amplitudes. In particular, all the conditions analyzed here that engage 
transfer, involve only large wavelengths and are therefore fairly accurate. 
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FIGURES 

FIG. 1. Poincare plots on the homogeneous manifold for r = 1.0: Ss = 0.0001 in (a) and 

^3 = 0.1 in (b) 

FIG. 2. Time series of V < > for r = 1.0, Vg^ = 1.0, N = 128, and various values of ^3. 
Initial conditions discussed in the text. 

FIG. 3. (a) Comparison of time series of V < N'^ > for different number of modes used in the 
simulations; (b) Spectral-energy distribution of the case N = 256 above, averaged over time and 
over the mode index j. In all cases 5s = 0.1, r = 1.0, Vg^ = 1.0; initial conditions discussed in the 
text. 

FIG. 4. Space-time history of |^i(a;,t)| for the 63 = 0.1 case of Fig. (2). By the end of the 
simulation inhomogeneities are already strongly excited. 

FIG. 5. Time series of V < TV^ > when orbits of the regular type evolve on the homogeneous 
manifold. The number of linearly unstable modes in all cases is approximately 5. In all cases 
S3 = 0, Vg2 = 0.06 and A'" = 128; initial conditions discussed in the text. 
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Fig.2: Lopes & Rizzato 
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Fig. 3a: Lopes & Rizzato 
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Fig.Sb: Lopes & Rizzato 
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